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Abstract. It is widely believed that multiwavelength afterglows of gamma-ray bursts 

■ (GRBs) originate from relativistic blast waves. We here show that in such blast waves, 
r ~ l ' a significant fraction of the energy of shock-accelerated protons would be lost due to 
I/"*) ' pion production by interactions with afterglow photons. This could lead to long-term 

m '. 

■^J- , production of 10 16 -10 18 eV neutrinos and sub-TeV 7-rays that accompany with usual 

■ afterglows, provided that the protons are accelerated to 10 19 eV in the blast waves. 
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1. Introduction 

In the standard fireball shock model, GRBs and their afterglows are produced due to the 
dissipation of kinetic energy of relativistic ejecta with Lorentz factor of 70 ~ 100 — 1000, 
so-called "fireball", released from an unknown central engine (see Piran 1999; Paradijs, 
Kouveliotou, & Wijers 2000; Meszaros 2002 for detailed reviews). In this model, the 
shocks, resulting from the collisions between shells with different Lorentz factors in the 
ejecta (internal shocks) or between the ejecta and the surrounding medium (external 
shocks), accelerate the electrons responsible for prompt gamma-rays in GRBs and for 
X-ray, optical and radio emissions of afterglows. The same shocks should also accelerate 
protons. It is pointed out that the physical conditions in the shocks allow protons to be 
accelerated to e p ~ 10 20 eV (Waxman 1995a; Vietri 1995; Dermer & Humi 2001. See 
but Stecker 2000; Mannheim 2000; Scully & Stecker 2002). Furthermore, the spectrum 
and flux of ultra-high energy cosmic rays (UHECRs) are consistent with those expected 
from Fermi acceleration of protons in GRBs (Waxman 1995b), suggesting a common 
origin of GRBs and UHECRs. This picture will lead obviously to high energy neutrinos 
produced by ir + created in photo-pion interactions between the accelerated protons in 
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the internal or external shocks and synchrotron photons from electrons accelerated by the 
same shocks (Waxman & Bahcall 1997, 2000; Halzen 1998; Vietri 1998ab; Dai & Lu 2001; 
see Waxman 2001 for a review). This pj reaction has a peak at the energy threshold of the 
photo-meson A-resonance, i.e. e 7 e p = 0.2GeV 2 7 2 , and the resulting neutrinos own ~ 5% 
of the proton energy in the decay of charged pions and muons. This leads to neutrinos 
of ~ 10 14 cV for prompt gamma-rays from internal shocks (Waxman & Bahcall 1997; 
Halzen 1998), while a UV/optical flash from the external reverse shock implies higher 
energy neutrinos (e.g., - 10 18 eV) (Waxman & Bahcall 2000). Vietri (1998a, 1998b) 
studied the neutrino production in the case of external forward shocks, and pointed out 
that the neutrino energy may exceed 10 19 eV. But Waxman & Bahcall (1999) argued 
that the contributions of external shocks to the neutrino flux is small. 

In recent years, the standard afterglow model, which is significantly successful in inter- 
preting observations, has been developed (Meszaros & Rees 1997; Sari, Piran & Narayan 
1998; Wijers & Galama 1999). In this paper, we re-investigate the case of external shock 
neutrinos based on the present detailed afterglow model. We derive the maximum energy 
that protons could reach by Fermi acceleration and the efficiency of energy loss from ac- 
celerated protons to pions in the forward shock. By considering the cooling of secondary 
particles, ir + and we derive the spectrum and the flux of the resulting neutrinos. 
We find that at late times, e.g., ~ 1 day after the GRB trigger, the forward shock may 
be able to produce an intense neutrino signal, even though the flux is smaller at earlier 
times, as stated by Waxman & Bahcall (1999). 

We discuss in section 2 the standard afterglow model, including the dynamics and 
synchrotron radiation of the blast wave, and in section 3 the proton acceleration by 
the same shock. In section 4 we discuss the neutrino production rate from the shock- 
accelerated protons and the maximum energy that neutrinos may achieve. The neutrino 
spectrum, flux and detectability are derived in section 5. Conclusions and discussion are 
given in section 6. 



2. Dynamics and emissions of afterglows 

The fireball shock model suggests that after the prompt burst from internal shocks, the 
ejecta with total kinetic energy E continues to expand and drives a relativistic blast wave 
into the surrounding medium. Once the shocked medium gets an energy similar to E, 
the blast wave will enter the Blandford-McKee (1976) self-similar evolution, where the 
radius r and the Lorentz factor 7 evolve with observer time as (Wijers & Galama 1999) 

r = 4.07 x 10 17 (S 53 /n ) 1/4 iy 4 cm, (1) 
7 = 8.87(i?53/no) 1/8 ^ 3/8 , (2) 



where n is the baryon number density of the medium, and td — t/1 day. Unless stated 
otherwise, the useful form Q = \Q X Q X has been used in this paper. 

Let £ e and £b be the fractions of the post-shock energy density U — 4^ 2 nm p c 2 (in the 
comoving frame with the shocked material) that are carried by electrons and magnetic 
fields respectively. The characteristic electron Lorentz factor (in the comoving frame of 
the shocked material) is j m ~ £,e"/rn p /m e , where the post-shock energy per proton is 
7m p c 2 . The strength of the magnetic field in the comoving frame is 

B = 0Me B %Eli 8 n 3 Q /s t- 3/8 G. (3) 

From the above equations, we can calculate the evolution of synchrotron spectrum 
from the blast wave. The characteristic energy of synchrotron photons from electrons 
with Lorentz factor 7 m is 

et - frr&jg = O.O^^^V^V, (4) 

where and hereafter we denote the particle energy in the observer frame with super- 
script "ob" , and particle energy measured at the comoving frame without superscript. 
The characteristic energy of synchrotron photons from electrons with the cooling time, 
6irm e c/ ot^IcB 2 , which equals the blast wave expansion time (~ r/^/c), is given by 

eB 

m e c 

and the specific luminosity L € = dL/de oh at the peak is 



t 2 e -° ot-3/2 -1 £,-1/2,-1/2 t t /r\ 
n lTc-—=^B,-2 n E 53 t d eV > ( 5 ) 



L max ~ (2 1 rn)~ 1 1 ^N e ~ 5 x 10 B8 (&,,_ 2 no) 1/2 .E53 s~\ (6) 
ui e c 

where N e = ^irr 3 n is the total number of the swept-up electrons. Eq. (^) implies that the 
peak specific luminosity is a time-independent constant. The electrons are accelerated by 
the shock to a power-law energy distribution, dN e /drf e oc j~ p for j e > j m , where p ~ 2.2 
is fitted to the observations of afterglows. We take p = 2 for simplicity here. After t ~ 1 
hour, e° b > e^ 5 . In this case, the specific luminosity L £ peaks at e*, L max = L m , 
and the synchrotron spectrum is a broken power law (Sari, Piran & Narayan 1998) 
as: L e — L m (e/e m )^ 1 ' 2 for e m < e < e c and L c oc e" 1 for e > e c . Note that in the 
above equations, the typical parameter values have been adopted as £ e = 10 £ ej _i, 
£b = 1CP 2 £b.-2, E = 10 53 £^53 erg and n = Ihq cm -3 , which arc inferred from some 
afterglow observations (e.g. Wijers & Galama 1999; Granot, Piran & Sari 1999). 

3. High-energy proton production 

Following the analysis of Waxman (1995a) in the case of internal shocks, we investigate 
the physical condition in external forward shocks for high energy protons production. The 
typical Fermi acceleration time is t a — fR^/c, where the Larmor radius is given by Rl = 
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Cp/eB and / is of order unity (Hillas 1984). The protons are available for acceleration 
only within a comoving expansion time td ~ r/'yc. This leads to the requirement that 
t a < t d . From Eqs. (0)-(||), we have 

£ ° b < 4 x lO 19 ^^ 7 V V /8 eV. (7) 

If a magnetic field is nearly in equipartition (i.e., £b ~ 1), the protons can be accelerated 
to well above 10 20 eV. However, fits of afterglow observations give a typical value £b = 
0.01, showing that the magnetic fields are too weak to accelerate protons up to 10 20 eV 
in the blast waves. This upper limit is somewhat robust since its dependence on the 
medium density and the time (oc n 

1/8^-1/8) 

is very weak. 

The proton acceleration is affected by energy loss due to synchrotron radiation. The 
requirement that the proton synchrotron cooling time, t PjSyn = {(Sirm^c^ / aTm~)€ p B~ 2 , 
should be larger than the acceleration time t a is 

e? < 3.2 x l^Cs^ET^ ^ V A6 eV. (8) 

The proton acceleration is also limited by energy loss due to interaction with afterglow 
photons. In the next section, we show that, provided a typical medium density uq = 1, 
the fraction of energy loss of a proton in the expansion time td is usually not larger than 
one, which implies that the limit from energy loss due to photo-pion interaction is a less 
stringent constraint than that from t a < td- Thus, the maximum energy of protons comes 
from Eq. ©, e° b mQ ^ 4 x 10 19 eV. 

Gallant & Achterberg (1999) studied the shock-acceleration, also called first order 
Fermi-acceleration, in GRBs and argued that it may be impossible for protons to be 
shock-accelerated to above 10 18 ' 5 eV in an ultra-relativistic blast wave which is driven by 
the fireball ejecta into a typical interstellar medium. However we still expect protons to 
be accelerated to ultra-high energy in the blast waves by second order Fermi-acceleration, 
where the protons gain energy from the scattering at plasma waves (e.g. Dermer & Humi 
2001). 

4. Neutrino production 

Denoting the photon number density in the comoving frame of the shocked medium by 
n 7 (e)de and following Waxman & Bahcall (1997), we can write the fractional energy- loss 
rate of a proton with energy e p due to pion production, 
1 de n 



= — c / deav(e)£(e)e / dxx~ 2 n(x), (9) 

where 7 p = t p /in p c 2 , e>v(e) is the cross section for pion production for a photon with 
energy e in the proton rest frame, £(e) is the average fraction of energy lost to the pion, 



eo = 0.15 GeV is the threshold energy, and the proton number density is related to the 
observed specific luminosity by n{x) — L e ( r yx) / (Airr 2 cjx) . Because of the A-resonance, 
the photo-meson production for eo/2e c *C 7 P < eo/2e m is dominated by interaction with 
photons in the energy range e m < e -C e c , where L e oc e -1 / 2 . For this photon spectrum, 
Eq. (0) leads to 

^_X/ \ ^ ^ ^ -^m Z' ^pcak \ ^"pcak^peak A6 (10) 

71 p 2.5 47rr 2 7 V7p £ ™/ e poak 

Here, <7p C ak — 5 x 10~ 28 cm 2 and £ pca k — 0.2 at the resonance e = e pca k = 0.3 GeV, and 
Ae ~ 0.2 GeV is the peak width. The fraction of energy loss of a proton with observed 
energy e° b by pion production, fir(tp h ), is estimated by t~ 1 (e p ) times the comoving ex- 
pansion time of the shocked material (~ r/^fc). Thus, according to Eqs. (Q), (Q), (Q) and 
(^), we obtain 

Uef) « om^e B %EU\l%^ 9 )^. (11) 

Eq. (|l"l| ) is valid for protons in the energy range 

6 x 10^l%EU 4 nrr d 1/4 eV < ef 
< 5 - ^X^lE^n- 1/4 ifeV . (12) 

At lower energy, protons interact mainly with photons in the energy range e > e c , where 
L e oc e _1 , and oc e° b . Eqs. (|ll]) and fll2] ) imply that the fraction of energy loss, /^h, at 
the "break" energy e° b = 6 x 10 17 -2-^53 n o *d 1 e ^ > * s a time-independent constant, 

/ w6 « 0.02e e ,-i^i 2 £53^o /2 - (13) 

The maximum energy of the resultant neutrinos is analyzed as follows. This energy 
is first determined by the maximum energy that the accelerated protons achieve. From 
Eq. (0), the maximum energy of neutrinos is ~ 2 x 10 18 eV. This energy is also limited 
by the energy loss of secondary particles, i.e., pions and muons, because both pions and 
muons may suffer synchrotron and inverse-Compton (IC) losses before decay (Rachen & 
Meszaros 1998). For synchrotron losses of pions and muons to be negligible, their Lorentz 
factor (in the comoving frame of shocked medium) should not exceed the critical Lorentz 
factor 7* given by 7 2 t* = r yptp^y n (m ir /m p ) 3 (Rachen & Meszaros 1998), where * refers to 
either pions or muons, t„ = 2.6 x 10~ 8 s and r M = 2.2 x 10 -6 s are respective lifetimes in 
their rest frames, and m^/rrip ~ 0.1 for both pions and muons. Due to the longer lifetime 
compared to pions, muons suffer greater energy loss, resulting in lower upper limit of 
neutrino energy. Using Eq. (|3|), we have (in the observer frame) the critical energy for 
muons 

in = 10 2 °^ 2 2 n- 1/2 eV, (14) 
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which is independent of the total fireball energy E and observer time t. For IC loss to 
be negligible, the critical Lorentz factor is given by 7^-r* = r yptp ta y n (m i , /m p ) 3 (Ub /U p h) , 
where Ub = £bU and U p h — £, e r]U are the energy densities of the magnetic field and 
afterglow photons, respectively, and 7/ is the fraction of the electron energy that is radiated 
away. Note that 77 = (tc/fm)' 2 ^' 2 for slow cooling electrons with e c > e m . We take 
p « 2.2. 77 is not sensitive to other parameters as 77 ~ f. (In fact, 77 — I shows the 
strongest IC loss and provides the strongest constraint). Thus, the critical Lorentz factor 
is7^«7^ n (W£e) 1/2 ,or 

e ° b /c = 4xl0 19 C / iV /2 eV, (15) 

which is independent of £g. The IC loss is somewhat stronger than the synchrotron 
loss provided £s/£ e < 1- But, the energy limit from both losses exceeds the maximum 
proton energy from Eq. (Q), so the neutrino production does not suffer the energy losses 
of secondary particles. The maximum neutrino energy is determined by the maximum 
proton energy, e° b ~ 2 x 10 18 eV. 



5. Neutrino spectrum, flux and detectability 

The photo-meson interactions include (1) production of ir mesons: pj — > p + ir° and 
737 — » n + 7r + , and (2) decay of 7r mesons: ir° — > 27 and 7r + — > /j + + i/p — » e + + v e + 
*V + v n- These processes produce neutrinos with energy ~ 5%e p (Waxman & Bahcall 
1997). Eq. ([l2] ) implies that for a fixed time, the spectrum of neutrinos below w 
3 x 10 16 ^g _ 2 £^53 n.Q 4 t j X ^ 4 eV is harder by one power of the energy than the proton 
spectrum, and by half a power of the energy at higher energy. Therefore, if the differential 
spectrum of accelerated protons is a power law form n(e p ) oc e p 2 (Blandford & Eichler 
1987), the differential neutrino spectrum is n(e„) oc e^ 1 below the break and n(e„) oc 
e„ above the break. Unlike the prompt bursts of neutrinos from either internal or 
reverse shocks, which has a short duration similar to the GRBs, the neutrino emission 
from forward shocks last for a longer time due to the long-term forward shocks. Without 
the detailed knowledge of the evolution of the spectrum of accelerated protons in the 
forward shocks, it is far from calculating clearly the neutrino fluxes. However, as e°^(oc 
t~ 1//4 ) and £p b mQa: (oc t -1 / 8 ) do not change much from one day to even weeks which is 
concerned about here, one can expect a time-integrated neutrino spectrum similar to the 
spectrum at t = 1 day. For simplicity, the particle energies, hereafter in this section, refer 
to the value when t = 1 day, e.g. tf b -> tf b {t d = 1), e° b -> ef b (t d = 1), etc. 

We derive the neutrino fluxes as follows. Defining by f v = 
J n^e^epf^dep/ J n(ep)e p de p the fraction of energy loss of all accelerated protons 
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to neutrinos, and considering the neutrino spectrum given in the above paragraph, we 
obtain, after some calculations, 

/ob / e obU/2 
j , \ c p,max/^pbJ , ^ ob 

LlL \tp,rnaxl ^pl ) 

Here f„b is a constant for a certain GRB, therefore this value is mainly bound to e p ; 
which is determined by the lower threshold of a detector. In this equation, e°j? should be 
replaced by e p i for e pl > e° b . Taking e pi = e°£, we have f n « 0.07. 

The present day muon neutrino energy flux due to GRBs is approximately given by 
J w 0.25(c/4:ir)f n EtH, where tn ~ 10 Gyr is the Hubble time, and E is the production 
rate of UHECRs in GRBs per unit volume. The factor 0.25 here accounts for that about 
one half of the proton energy is lost to neutral pions which do not produce neutrinos, 
and about one half of the energy in charged pions is converted to v^s and v^s. We take 
E ~ 10 44 i?44ergMpc _3 yr _1 (Waxman 1995b). Integrating the specific neutrino flux 
over the neutrino spectrum (<&„ cx [e° b ]° below e°\, and oc [e° b ] -1 ^ 2 above e° b ), we 
have J = J <& v de v ~ 2(e° b e° b ma;r ) 1 / 2 $^t. Combining these two equations about J one 
obtains that the neutrino flux at is 

^«, xl(r »A^ _^-jS~- ) " 2 cm-s-.3 t -. ,17, 

The resulting high-energy neutrinos may be observed by detecting the Cherenkov light 
emitted by upward moving muons produced by neutrino interactions below a detec- 
tor on the surface of the Earth (Gaisser, Halzcn & Stanev 1995; Gandhi et al. 1998). 
Planned 1 km 3 detectors of high energy neutrinos include ICECUBE, ANTARES, 
NESTOR (Halzen 1999) and NuBE (Roy, Crawford & Trattner 1999). The probabil- 
ity that a neutrino could produce a high-energy muon in the detector is approximated by 
P v ^n ~ 10 _3 (e y /10 16 eV) 1 / 2 . We obtain the observed neutrino event rate in a detector, 

iVcvents(> e°b) = 27r/ ^P u ^def « 27r$, b P^ M (e° b )ln 

/ ob \"V2 

« 0.04^4 (^^rrPtv) km-V" 1 . (18) 

6. Conclusions and discussion 

After a prompt 7-rays burst, the fireball ejecta continues to drive a relativistic blast wave 
into the surrounding medium. The shock wave accelerates the electrons to give rise to 
observed afterglow by synchrotron radiation. The protons in the shock are also expected 
to be Fermi-accelerated at the same time. However, the maximum proton energy may not 
well exceed ~ 10 20 eV based on the typical parameters of afterglows. The interaction of 
the protons with synchrotron photons leads to a significant fraction of proton energy loss 
to neutrino production at energy 10 16 -10 18 eV. Though the synchrotron photon energy 
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density in the blast wave is much lower than in the reverse shock, the efficiency of proton 
energy loss in the blast wave is still comparable to that in the reverse shock, because of 
the long expansion time of the blast wave available for energy loss of protons. Since both 
shocks carry comparable energies (Sari & Piran 1999), we assume the same production 
rate of high-energy cosmic rays, E44 ~ 1 for blast waves as Waxman & Bahcall (2000) for 
reverse shocks. Thus, the expected detection rate of muon neutrinos 0.04 km~ 2 yr^ 1 is 
also comparable to the reverse shock. These mean that besides the internal shock and the 
reverse shock, the blast wave is also one region suitable for intense neutrino production. 

Waxman & Bahcall (1999) have estimated a much lower efficiency of the blast wave, 
i.e. ~ 10~ 4 , using Eq. (4) of Waxman & Bahcall (1997). However, the parameter 
values they chose are only valid for the early phase when the external shock just forms. 
When a set of parameters, which are valid for afterglows one day after the burst, are 
chosen as follows: L 1 = 10 45 erg/s, At = 10 5 s, e° b = O.leV and T = 10, their equation 
gives an efficiency f v ~ 10~ 2 in general agreement with ours by orders of magnitude. 

It may be that UHECRs can not be produced in an ultra-relativistic forward shock 
(Gallant & Achterberg 1999) during the transition of the shock to the self-similar ex- 
pansion. However the hard X-rays from the forward shock will interact with the protons 
accelerated by the reverse shock, yielding a new component of neutrino emission besides 
those produced by interaction between protons and optical/UV photons from the reverse 
shock, which is proposed by Waxman & Bahcall (2000). At the moment of transition, 
especially for the photons at energy e ob <; ~ 80 keV (cf. Eq. (J]), the specific lumi- 
nosity from the forward shock dominates that from the reverse shock, which implies a 
higher neutrino flux at energy e° b £ 10 15 (7 /250)(e° b /100keV)- 1 eV. 

Different models for progenitors of GRBs include different environments: compact 
object mergers occur in the ISM and the explosions of massive stars in the preexisting 
stellar wind or denser medium. Dai & Lu (2001) have found that the neutrino production 
from reverse shocks in the wind case has a quite different spectrum and 10 times more 
detection rate compared with the ISM case, which provides a way to distinguish between 
GRB progenitor models. In the blast-wave case discussed in this paper, the proton energy 
loss rate shows a strong dependence on density n(Eq. we expect much higher fluxes 

for a denser medium. For 714 — 1 typically for giant molecular clouds, following the 
analysis in this paper, the energy loss efficiency of a proton at one day after the the GRB 
trigger is given by f„ w 8.6e° b 18 ty 4 for e° b < 5 x 10 18 i^ 4 eV, implying that protons of 
energy e° b > 10 17 eV lose all their energy to pion production and the neutrino spectrum 
for e° b > 5 x 10 15 eV will trace the proton one. The maximum energy that protons can 
be accelerated to, which is now confined by the pion production, is 4 x 10 18 ^ 3/4 eV. 
This implies a maximum neutrino energy of ~ 10 17 eV. 
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Another consequence of the photo-meson interactions is a sub-TeV afterglow. The 
decay of neutral pions produced in the blast waves would lead to ~ 10 18 eV photons. Such 
high energy photons would suffer pair production in the afterglows because the optical 
depth for photons above ~ 1 TeV is important (Zhang & Meszaros 2001). So these ultra- 
high-energy photons would be degraded and leak out as a long-term sub- TeV afterglow, 
carrying away ~ 10% energy of protons, similar to the energy losses to neutrinos. For 
a GRB at typical redshift z = 1, the intergalactic diffuse radiation field absorbs the 
photon above 100 GeV (Mannheim, Hartmann & Funk 1996), allowing only a multi-GeV 
afterglow observed on the earth. An EGRET observation of GeV emission extended to 
1.5 hours after the trigger of GRB 940217 (Hurley et al. 1994) has indicated that a high 
energy spectral component can extend to such high energy band and persist for a long 
period of time. However, the electron-IC emission in the blast wave can also contribute to 
this GeV observation. For the parameters, no — £e,-i — £,b,-2 — 1, adopted in this paper, 
the electron-IC emission is important and plausible to explain the observation (Zhang & 
Meszaros 2001). From the afterglow model presented in section 2, we can see that the 
peak energy in the IC emission spectrum, 7^e° b (Sari & Esin 2001), is still marginally up 
to 100 GeV even one day after the GRB trigger, but we do not expect GeV IC emission a 
few days later, while the emission produced by neutral pion decay might still be possible. 
In principle, the GeV emission could also be due to the proton synchrotron emission, but 
it is not important for the parameters adopted here (Zhang & Meszaros 2001). 
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